The curvaton scenario within the MSSM and predictions for non-Gaussianity 
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We provide a model in which both the inflaton and the curvaton are obtained from within the 
minimal supersymmetric Standard Model, with known gauge and Yukawa interactions. Since now 
both the inflaton and curvaton fields are successfully embedded within the same sector, their decay 
products thermalize very quickly before the electroweak scale. This results in two important features 
of the model: firstly, there will be no residual isocurvature perturbations, and secondly, observable 
non-Gaussianities can be generated with the non-Gaussianity parameter /nl ~ 0(5 — fOOO) being 
determined solely by the combination of weak-scale physics and the Standard Model Yukawas. 



The curvaton scenario [1-4] is an alternative 
mechanism for the generation of the primordial per- 
turbations whose spectrum is observed in the cosmic 
microwave background (CMB) Q. In this scenario, 
the density perturbations are sourced by the quan- 
tum fluctuations of a light scalar field 0, the cur- 
vaton, which makes a negligible contribution to the 
energy density during inflation and decays after the 
decay of the inflaton field a. (For a review on infla- 
tion including the curvaton mechanism, see [6] .) The 
advantage of the curvaton mechanism is that it can 
in principle generate measurable non-Gaussianity 
[l|,|7|in the primordial density perturbations and also 
significant residual isocurvature perturbations, nei- 
ther of which are possible in the usual single-field in- 
flation models. Both signatures are detectable, and 
if either were to be observed, this would strongly 
favour the curvaton hypothesis. 

If the curvaton does not completely dominate the 
energy density at the time of its decay, the process of 
conversion of initial isocurvature perturbations into 
adiabatic curvature perturbations can enhance non- 
Gaussian fluctuations to the level where they might 
be constrained by the Planck satellite. The enhance- 
ment in non-Gaussianity is given by /nl ~ 5/(4r) 
for r < 1, where r = Pcf,/ p-cg.^ at the time the cur- 
vaton decays [l|. Planck is expected to be able to 
detect non-Gaussianity of the order /nl ^5 [8]. To 
achieve detectable /nl thus requires small r. 

However, if either the curvaton or the inflaton be- 
long to a hidden sector beyond the Standard Model 
(SM), they may decay into other fields beyond the 
SM degrees of freedom [dof). There is no guaran- 
tee that the hidden and visible sector dof should 
reach thermal equilibrium before Big Bang Nucle- 
osynthesis (BBN) takes place. In this case, resid- 
ual isocurvature perturbations are expected to be in 
conflict with CMB data, which constrain them to be 
less than 10% [^. If the curvaton belongs to the vis- 
ible sector but the inflaton does not, a value of r '^ 1 
would avoid this conflict [10[ but would render any 
non-Gaussianity undetectable. Note that if r ^ 1 



the curvaton is solely responsible for exci ting all the 
SM dof so it must carry the SM charges lll.[l2|. 

For the curvaton model to be observationally dis- 
tinguishable, we wish the model to be able to cre- 
ate detectable non-Gaussianity. For this, r must 
be small and both the inflaton and curvaton decay 
products must thermalize before the time of nucle- 
osynthesis, as there are stringent constraints on any 
non-SM like hidden radiation after BBN [Oj. In order 
to achieve this, we wish to place the entire inflaton- 
curvaton paradigm within a particle physics model 
where all the interactions are well constrained by the 
weak scale physics. 

Recently, the inflationary paradigm has been em- 
bedded within the minimal supersymmetric Stan- 
dard Model (MSSM) with known gauge interac- 
tions [ij, [ij] . The aim of this letter is to show, for 
the first time, that it is possible to embed both the 
inflaton and curvaton within MSSM, without involv- 
ing any hidden sector. We thus provide a solution 
to a general problem of the curvaton scenario, i.e., 
how to generate measurable non-Gaussianity with- 
out large residual isocurvature fluctuations. 

Let us flrst consider the total potential to be the 
sum of inflaton vacuum energy, denoted by Vb , and 
curvaton potential V{(j)) 



Vtotal ^Vo + V{(j)) . (1) 

We assume V"{(l)) ~ r7i^((/'/) < i/| - Vo/M^ 
(Mp ^ 10^® GeV) where the subscript / indicates 
the quantities are evaluated during inflation. This 
condition is required for a successful curvaton sce- 
nario. The curvaton acquires vacuum induced quan- 
tum fluctuations, which have amplitude 
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These fluctuations are converted into the adiabatic 
density perturbations when the curvaton decays dur- 
ing its coherent oscillations or rotations. In order to 
match the observed amplitude of the fluctuations on 
the GMB, r6 - 10"^ 



Let us first discuss the origin of tlie curvaton, 
wfiicfi we take to be an i?-parity conserving D-fiat 
direction of the MSSM (for a review see [lH). Two 
candidate flat directions are LLe (where L denotes 
the left-handed slepton superfield and e the right- 
handed superfield) and udd (where u and d de- 
note the right-handed squark superfields), which are 
lifted by the non-renormalizable operator: 
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where A is a non-renormalizable coupling. For con- 
creteness, we take the curvaton to be LLe so that 
the scalar component of the $ superfield is: 



= (i + L + ?)/%/3, 



(4) 



where L and e are the slepton and selectron scalar 
fields. At the lowest order the potential along the </> 
direction is given by: 
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(5) 

where A - m^ ~ 0(100-1000) GeV, m^ is the soft 
SUSY-breaking mass term, and n — 6 for LLe [15| . 



During inflation if mi <C 
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the fluctuations 



along this nearly massless direction would create a 
homogeneous condensate with a vacuum expectation 
value (VEV) given by [H] 



(j)i - {m^M^- 
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10^* GeV 



(6) 



assuming A - 0(1). For m^ - 100 - 1000 GeV, 
and n — 6, in order to match the amplitude of the 
density perturbations 6, the Hubble expansion rate 
during inflation should be Hj ^ 10^° GeV if r ^ 1. 

There is a distinction between a positive and neg- 
ative phase of the A term. The difference in dynam- 
ics arises after the end of inflation. In the case of 
positive A term the curvaton starts rolling towards 
the origin immediately, but in the case of a negative 
phase, for values of ^ > v4Or?T,0, it may remain in a 
false vacuum with the VEV given by Eq. (j6|). In this 
case the curvaton rotates instead of oscillates around 
its global minimum at cf) — 0. In either scenario, the 
curvaton mass is negligible compared to the Hubble 
expansion rate. In fact, for A = ^/40m^ and a neg- 
ative phase the curvaton is actually massless along 
the real direction, and obtains inflaton-induced ran- 
dom fluctuations of order Scj) ~ Hi/2tt. 

We now turn to the origin of Vq within the MSSM. 
Let us consider a flat-direction orthogonal to the cur- 
vaton. If the curvaton is LLe, this could be the udd 
direction. We take the inflaton direction to be: 



a^ {u + d + d)/Vi, 



(7) 



where u and d are squark scalars. Note that udd 
and LLe remain two independent directions for the 
entire range of VEVs. 

This flat direction will also be lifted by the non- 
renormalizable operators. However, at larger VEVs 
the potential energy density stored in the udd di- 
rection will be larger than for the LLe, so it would 
be lifted by higher order terms: 
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The potential at lowest order would be: 
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where . . . contain the higher order terms. Note that 
the Am in Eq. ([5]) are all non-renormalizable cou- 
plings induced by either gravity or by integrating 
out the heavy fields at the intermediate scale. At en- 
ergies below the cut-off scale these coefficients need 
not necessarily be of 0(1). 

Potentials like Eq. ([SJ were studied in Refs. [la, 
[13. For A2 < A3 <C A4 < A„ < 0(1), they 
provide a unique solution for which first and sec- 
ond derivatives of the potential vanish along both 
radial and angular direction in the complex plane: 
dV/da = dV/da* = d^V/da^ = d^V/da*^ == (a 
saddle point condition) [18] . For the first three terms 
in Eq. ([9]), it is possible to show that this happens 
when 
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at the VEVs: a = ctq exp [?7r/3, in, i57r/3], ctq = 
(2/ll)(A3/A4)^/^Mp. Concentrating on the real di- 
rection, the potential energy density stored in the 
inflaton sector is given by: 
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(11) 



where (To ^ Mp . Note that inflation occurs near the 
saddle point ao, where the effective mass vanishes. 
However, the third derivative of the potential is not 
negligible, V" -- Xjcrl/M^ ^ 0, which leads to slow 
roll inflation. The potential for the inflaton becomes 
flat enough to sustain a large number of e-foldings. 

As written, the condition Eq. (|10l) represents a 
complete fine-tuning. Some deviation from this con- 
dition will be possible, changing the saddle point to 
a point of inflection, so long as V remains small 
enough for sufficient e-foldings of inflation. Detailed 
discussion on fine-tuning in infiection point inflation 
can be found in Rcf. [13j |. 

The amplitude of perturbations of the inflaton 
is given by Sh ^ V"'icro)N^oBE/^07rHi ^. The 



corresponding Hubble expansion rate is given by 
Hi - (153/88)A2cr^/M^. For ctq ~ lO^^'^ GeV and 
A2 ~ 10~6, it is possible to obtain Hi - 10^° GeV, 
required for a successful curvaton scenario. For the 
above values, the inflaton perturbations are negligi- 
ble, i.e. 5h < 10~^, therefore all the observed per- 
turbations are created mainly by the decay of the 
curvaton. 

Now let us consider the aftermath of inflation. 
The inflaton would decay primarily into the MSSM 
dof. The coherent oscillations of the inflaton would 
give rise to instant preheating and thermalization of 
the light MSSM dof as discussed in Ref. [19], with a 
reheat temperature 



Tr ^ [HiMp] 



1/2 



10" GeV . 



(12) 



However, not all of the MSSM dof will be in ther- 
mal equilibrium in our case. For the given choice 
of flat-direction fields, if both inflaton and curvaton 
simultaneously take large VEVs, the SU{2)w dof 
would not reach in thermal equilibrium, since the 
LLe VEV would induce large masses to those dof. 
This will play a crucial role in determining the non- 
Gaussianity parameter /nl, as we shall show below. 

The curvaton (j) starts to rotate about the ori- 
gin when H = Hose "^ m^. The field value at 
this time is |0osc| ^ (m^Mp^^y/"'-^. During this 
epoch the universe is already radiation-dominated 
following the decay of the inflaton. However, the 
curvaton cannot decay immediately, due to the fact 
that the curvaton VEV induces large masses h{<j){t)) 
for gauge bosons, gauginos and (s)leptons, where 
h is the gauge or Yukawa coupling. The curva- 
ton's decay at leading order is kinematically forbid- 
den if h{<j)) > m^/2 - 0(100 - 1000) GeV. Decays 
do not occur until the Hubble expansion has red- 
shifted {(j){t)) down to 7710/2/1. Note that the SM 
Yukawa couplings are smaller than the gauge cou- 
plings. Therefore the decays via SM Yukawas be- 
come kinematically allowed at higher VEVs. 

During the rotations, the curvaton VEV will scale 
as (/)(i) oc a~^/^, as a oc H~^/^ during the radiation- 
dominated epoch. Therefore, each decay channel be- 
comes allowed when [20| 

H = iJdcc - m^ {m^lh(j){t)f^ , (13) 

For large {4>{t)), the decay time is naturally longer 
than the normal decay rate into the massless dof. 
The radiation energy density stored in the infiaton 
decay products scales as pyis oc H^, where the sub- 
script denotes the visible dof. The ratio of the energy 
densities at the time the curvaton decays is given by 
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The kinematical blocking due to the curvaton VEV 
enhances the efficiency factor, r, therefore the cur- 
vaton rotations prolong the mater-dominated epoch 
till it decays completely. For soft SUSY-breaking 
mass rrirf, < 1 TeV, the inefficiency parameter is 
r ~ 0(l)h?/'^. Although the LHC has already placed 
severe constraints on the parameter space for low- 
scale SUSY, the current limits do not exclude heavy 
squark and slepton masses > 500 GeV [22|. Since 
our flat directions are all made up of squarks and 
slcptons, there is a large parameter space available 
in which this condition may be satisfied if SUSY is 
discovered at the LHC. 

Since the curvaton decay is delayed due to the 
kinematical blocking, r < 1 is different for each de- 
cay channel. What range of /nl we expect from the 
various dominant decay channels of the curvaton de- 
pends on the different values of h. If we consider the 
SM gauge couplings, then /i ~ 0.1 and we would ex- 
pect the largest /nl ~ (5/4r) - C'(l)/i-2/3 _ ^(5)^ 
However, the curvaton also has the Yukawa interac- 
tions, especially when the curvaton decays into lep- 
tons and sleptons, for which: 



/nl-^-O(I)/^ 



-2/3 



10- 10-^ 



(15) 
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for h ~ 10~^ — 10~^. This range of h covers all the 
SM Yukawas except the top Yukawa which is of order 
h ^ 0.1. Due to the smaller values of h, these decays 
are kinematically allowed at higher VEVs. An exact 
prediction for net effect on /nl requires a complete 
analysis of the decay modes for the LLe curvaton 
which is beyond the scope of the current letter, but 
it can be seen that this model of the curvaton can 
provide /nl in a range which will be observationally 
relevant in the near future. 

The temperature at which the curvaton decay 
products reach thermal equilibrium is determined by 
Eq. (Hg. The final thermal bath filled with MSSM 
dof would be obtained by the reheat temperature 

Trj ^ (ildecMp)'/' - 10*-' - lO^-^ GeV (16) 

for h ^ 10^^ — 10^'°^. Such a temperature is suffi- 
cient to excite weakly interacting massive particles 
and for baryogenesis |21j]. Note that both the tem- 
peratures from Eqs. (fT2|) and p6)) are sufficiently 
high to excite thermal/non-thermal gravitinos and 
axinos. If the gravitinos or axinos are the lightest 
SUSY particle, this causes two problems for this 
scenario: over-production of gravitinos with both 
helicities would be bad for BBN, and the graviti- 
nos and axinos would thermally decouple even be- 
fore the curvaton has started decaying. This would 
generate large residual isocurvature perturbations, 
because gravitinos and axinos can never come into 
thermal equilibrium. Instead the ideal dark matter 
candidate would be the neutralino, which decouples 
from the thermal plasma at T ^ 40 — 50 GeV. 



Our discussion so far has been based on treating 
udd as the inflaton and LLe as the curvaton flat 
direction. In principle, we could have swapped the 
roles of inflaton and curvaton, i.e. LLe as an infla- 
ton and udd to be the curvaton. The main aspects 
of the analysis would not differ at all. Although 
treating udd as a curvaton would also make SU{3)c 
dof heavy during the curvaton oscillations and this 
would alter the detailed discussion of thermalization, 
nevertheless the range of /nl quoted above for the 
SM Yukawas in Eq. (|15l) would remain the same. 

To summarize, we have discussed the possibility 
of constructing a model in which both the infla- 
ton and curvaton are flat direction flelds within the 
MSSM. The radiation created from the decay of the 



inflaton and curvaton belongs to the visible sector, 
avoiding the problem of residual isocurvature fluctu- 
ations, while the curvaton mechanism can create ob- 
servable non-Gaussianity. The non-Gaussianity pa- 
rameter /nl depends crucially on the SM gauge and 
Yukawa couplings, and ranges from 0(5) to 0(1000) 
in the different decay channels (for Yukawas in the 
range h ~ 10~^ — 10~^, which is the case for all 
the Yukawas except the top). The model favours 
a visible-sector dark matter candidate such as the 
lightest neutralino but will not work if the lightest 
SUSY particle is a gravitino or axino type. 
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